Carbon nanotube-based quantum pump in the presence of superconducting lead 



Yadong Wei^ and Jian Wang^'") 
1. Department of Physics, School of Science, Shenzhen University, Shenzhen 518060, China 
2. Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong, China 



O 
O 



0\ 



> 

m 
o 

O 

a 

o 
o 



Parametric electron pump through superconductor- 
carbon-nanotube based molecular devices was investigated. 
It is found that a dc current, which is assisted by resonant 
Andreev reflection, can be pumped out from such molecular 
device by a cyclic variation of two gate voltages near the nan- 
otube. The pumped current can be either positive or negative 
under different system parameters. Due to the Andreev re- 
flection, the pumped current has the double peak structure 
around the resonant point. The ratio of pumped current of N- 
SWNT-S system to that of N-SWNT-N system {I^^/l'^) is 
found to approach four in the weak pumping regime near the 
resonance when there is exactly one resonant level at Fermi 
energy inside the energy gap. Numerical results confirm that 
in the weak pumping regime the pumped current is propor- 
tional to the square of the pumping amplitude Vp, but in 
the strong pumping regime the pumped current has the lin- 
ear relation with Vp. Our numerical results also predict that 
pumped current can be obtained more easily by using zigzag 
tube than by using armchair tube. 

73.63.Fg,85.35.Kt,73.40.Gk,74.50.+r 



I. INTRODUCTION 

Physics of parametric electron pump has attracted 
great attention. Classical, pumps had been fabri- 
cated about a decade agaJ'cl. Recently, Quantum-dot- 
based quantum pump has bfijep- the subject of both 
experimentaH and theoreticalo c3 investigations. The 
quantum pump generates current due to the cyclic varia- 
tion of at least two system parameters while maintaining 
zero bias. As the charge pumped out of the system, it 
also produces the Joule heat along with the dissipation. 
Recently, Avron et alEil have derived the lower bound for 
the dissipation. This naturally leads to the concept of 
optimal pump. To search for an optimal pump, the heat 
current and the shot noise of quantum pump have been 
investijgated using the time-dependent scattering matrix 
theorynjcj. Veryjiecently, to explain the experimentally 
observed anomalyrl,|_th£L finite frequency pumping theory 
has been developedoE^. The adiabatic pumping theory 
has also been extended to account for the Andxe^v reflec- 



tion in the presence of superconducting 



I, strong 



electron interaction in the Kondo regimecS, and spin po- 
larized |pumped current when the ferromagnetic leads are 
presentEj. Due to the pcpuliar electronic properties of 
Carbon nanotube (CNT)Era, CNT-based parametric 
electron pump has been investigated as a prototypical 
nanometer-scale molecular devicetil. It would be interest- 
ing to further explore how does Andreev reflection modify 
the quantum interference of CNT based quantum pump 
in the presence of superconducting lead. It is well known 
that in the presence of normal conductor-superconductor 
(NS) interface, an incoming electron-like excitation can 
be Andreev reflected as a hole-like excitationEj. In this 
paper, we will study a hybrid structure where both car- 
bon nanotube and superconducting lead are present and 
examine the interplay between the electronic properties 
of CNT and superconductivity. Specifically, we investi- 
gate a parametric quantum pump that consists of a finite 
sized single wall carbon nanotube (SWNT) connected to 
one normal left lead (N) and one superconducting right 
lead (S), i.e., the parametric quantum pump of N-SWNT- 
S system. Two pumping driving forces Xi(t) and X2{t) 
are established by applying cyclic, time-dependent volt- 
ages to two metallic gates, which are capacitively coupled 
to the SWNT (see inset of Fig.l). We teve used nonequi- 
librium Green's function approachcB^Ea to calculate the 
pumped current. We found that in the presence of su- 
perconducting lead the pumped current is four times as 
that of corresponding normal system in the weak pump- 
ing regime. As the pumping amplitude increases, the de- 
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pendence of pumped current crosses over from quadratic 
to the linear dependence. Due to the Andreev reflection, 
the pumped current exhibits double peak structure for 
single resonant level in line with the chemical potential 
fig of the superconducting lead. When /i^ is in the middle 
of two resonant levels, another type of Andreev reflection 
occurs where an electron coming from normal lead tun- 
nels via the lower resonant level and Andreev reflected 
as a hole through the upper resonant level with a Cooper 
pair created in the superconducting lead. For this two 
level Andreev reflection, the pumped current shows re- 
markable parity effect that the direction of pumped cur- 
rent near one resonant level is opposite to that of the 
other level. Finally, our numerical results show that it 
is much easier to pump current through zigzag structure 
than armchair structure. 



II. THEORY 

In this work, we assume that the variation of the pump- 
ing potentials are very slow and the adiabatic approxima- 
tion is appropriate. In this approximation, the pumped 
current flowing through the left.. normal metallic lead, in 
one cycle time r, is given byQE3c3 
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where the quantity dN^/dXj is th£j)artial density of 
states (DOS), called the injectivitycjifj, of the left lead 
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where j — 1 or 2. dN^/dNj describes the number of 
electrons coming from left lead and exiting the system as 
electrons due to the external parameter Xj. dN^/dNj 
describes the number of holes coming from left lead and 
exiting thepSystem as electrons due to the external pa- 
rameter ATjO. G|Hi and 6*12 are the matrix elements of 
t^^-|2 X 2 Nambuo representation and can be expressed 
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Once the electron and hole Green's function G^i and 
were obtained, G12 is calculated by 
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is the line-width function and 



is the total self energy given by 
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where Y7^ = Pa — iTal2 is the self energy of the lead a 
in the normal case. Here Pa is the real part and Fq, is 
the linewidth function. The self energy for the supercon- 
ducting lead is 



Pr - iTRf3il2 
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where /3i = vE/^/E"^ - A^, = vAj^jE'^ - with 
1/ = 1 for > —A and = — 1 for £' < —A. Here 
A is the gap energy of superconducting lead and chem- 
ical potential of superconducting lead has been set 
to zero. In the above equations, is the Hamilto- 
nian of CNT. It is a x matrix, where A^ is the 
total number of carbon atoms. Vpp is a diagonal ma- 
trix describing the variation of the CNT potential land- 
scape due to the external pumping potentials Xi(t) and 
X2{t). In this work, we choose the two pumping po- 
tentials to be Xi{t) = Vi{t) = Vio + VipSin{u;t) and 
X2{t) = V2{t) = V20 + V2p sin{ujt + (p), where is phase 
difference, uj is the pumping frequency and Vip and V2p 
are the pumping amplitudes. To simplify the numerical 
calculation, we mimic the gate effects by simply adding 
the potential Vpp = ViAi + V2A2 to the SWNT where 
Ai is the potential profile function and can be set to be 
unit for the gate region, zero otherwise. A more accurate 
study requires a numerical solution of the Poisson equa- 
tion with the gates providing the appropriate boundary 
conditions. 



III. RESULT AND DISCUSSION 

We now apply Eq.(l) to calculate the current for the 
N-SWNT-S parametric pump. For simplicity, the SWNT 
is modeled with the nearest neighbor 7r-orbital tight- 



binding model with bond potential Ypp^^ = — 2.75eV. 
This model gives a reasonable, qualitative description 
of the elect.mnic and transport properties of carbon 
nanotubesOCJ. RecenlJji, a S-SWNT-S device has been 
studied experimentally^. By tuning the transparency 
of the device, clear signals of Andreev reflection were 
observed. The dependence of the Andreev current on 
the device transparency, the behavior of the differen- 
tial resistance in the sub-gap region, as wejU, as the ob- 
served low-temperature resistance anomalyEil can all be 
explained theoretically using the 7r-orbital tight-binding 
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modeB. We assume that the SWNT is weakly coupled 
to the electrodes so that the pumping process is medi- 
ated by the resonant tunneling. We further assume that 
strong electron-electron interactions may be neglected. 
Without losing generality, we set the energy gap of the 
superconducting lead to be lAbmeV (the gap of Nio- 
bium). We_also apply the wide bandwidth limit for the 
self-energycj and consider the symmetric pumping, i.e., 
Vio = V20 = Vq and Vip — V2p = Vp. In the absence 
of pumping we have Vp = which forms a symmetric 
double barrier with barrier height Vq in the finite size 
nanotube. As a result, the discrete resonant levels are 
established within the double barrier structure. By ad- 
justing Vq, we can control the positions of energy levels 
inside the energy gap A. Since the pumped current is 
proportional to uj, we set uj = I for convenience. We 
also set h = 2m = e = 1. When pumping frequency 
UJ = lOOMHz, which is close to the frequency used in 
Ref. ^, the unit for the pumped current is 1.6 x lO^^^A. 
Finally, we do not consider the finite temperature effect 
and hence set temperature to zero. 

First, we consider an armchair (5,5) SWNT with 93 
layers of carbon atoms (total 930 atoms) . The two metal- 
lic gates that provide the pumping driving forces are lo- 
cated near the two ends of the SWNT from 10th to 28th 
layer, and 66th to 84th layer. We have chosen Vq « 2.75V 
so that there is only one resonant level in the energy 
gap and the level is aligned with the chemical potential 
of the superconducting lead in the absence of pumping 
voltage 1^ = 0. Fig.l shows the transmission coefficient 
versus Fermi energy Ep. For comparison, the transmis- 
sion coefficient for corresponding normal system (when 
A = 0) is also plotted. We see that for normal system, 
the transmission coefficient (dashed line) has Lorentzian 
lineshape. In the presence of superconducting lead, the 
transmission peak (solid line) is flattened and narrowed. 
Fig. 2 depicts the pumped current I^^ versus the Fermi 
energy for different pumping amplitudes Vp with = 7r/2 
and Tl ^ Tji = 0.0136eT^. Several interesting observa- 
tions are in order. (1). the pumped current has large 
amplitude only near the resonant level showing clearly a 
resonant assisting behavior. This is because the pumped 
current is related to the global density of states (DOS) 
and the DOS of the system reaches its maximum near 
the resonant levelcl. (2). the amplitude of pumped cur- 
rent has double-peak structure around the resonant level. 
To understand this behavior, we plot the Andreev reflec- 
tion coefficient versus Ep at different times t during the 
pumping cycle in the right inset of Fig. 2. The Andreev 
reflection coefficient gives large value only around the 
two pumping instants: t — Sn/A and t = In/A because at 
such moments the energy level of the SWNT is just in line 
with the chemical potential of the superconducting lead 
so that an excitation of hole can be reflected when there 
is an incident election near the Fermi surface, and vice 
verse. At other moments, Andreev reflection coefficient 
is very small and contributes little to the current integral 
in the time cycle. In the inset, we just plot Andreev coef- 



ficients at certain moments: to = 37r/4, to = 0.95 x 37r/4 
(a little smaller than to) and to — 1.05 x 37r/4 (a little 
lager than to). We see that all three curves have double- 
peak feature (for the solid line the double-peak is barely 
seen). To understand this, we examine the Andreev re- 
flection coefficient which takes on the Breit-Wigner form 
near the resonant level _EoE3: 



TAiE) 



^1^2 



4[E^-E^, + ^y 



El[T- 



5Tf 
(10) 



where 5V = Vi - Fa and F = Fi Fj. Here Fi and 
F2 are linewidth functions which characterize the ability 
of electron tunneling through the left and right barri- 
ers, respectively. At the moment to = 37r/4, the barrier 
heights for the left and right barriers are Vq -I- {'\/2/2)Vp 
and Vo — {■\/2/2)Vp, respectively. In this case, we found 
that the position of the resonant level remains the same 
as when Vp = 0, i.e., the resonant level is still at Eq = 0. 
From Eq.(nO), we have 
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Eq.(|rT]) indicates that if Fi = F2 so that ST = 0, 
then resonant Andreev reflection occurs at £^ = with 
Ta(S = 0) = 1. On the other hand, if Fi > F2 so 
that (5F > 0, Ta takes on a maximum value at _E = 
but this maximum value is less than 1. Furthermore, if 
Fi < F2 such that ST < 0, Ta is characterized by two res- 
onant peaks with Ta = 1 at energies E± — ±\/ — F(5F/2. 
Since the left barrier is higher than the right one, we have 
Fi < F2. This explains why we have double-peak feature 
for the solid curve in the right inset of Fig. 2. At the mo- 
ment to = 0.95 X 37r/4 or to = 1.05 x 37r/4, the resonant 
level shift so that Eo is no longer zero. The double-peak 
structure can still be explained using Eq. (|l0|) . Similar 
discussions apply to the moment t = 77r/4. From the 
above discussion, it is clear that the double-peak feature 
of Andreev coefficient around to = Stt/A and to = 7tt/A 
is responsible for the behavior of the pumped current. 
(3). the pumped current increases as the pumping am- 
plitude Vp increases (compare solid line, dotted line, and 
dashed line in Fig. 2). (4). For comparison, we also 
plot the pumped current for N-SWNT-N system with 
Vp = 0.0014F (long dashed line in Fig.2). We see that it 
gives only one broad peak with small magnitude rather 
than two peaks. This is because for normal system the 
transmission coefficient does not give two peaks at any 
pumping moment and the double-peak feature for trans- 
mission coefficient is the intrinsic feature solely due to 
the Andreev reflection at the NS interface. Due to the 
quantum interference between the direct reflection (with 
amplitude Ed) and the multiple Andreev reflection (am- 
plitude Ea), the pumped current I^^ of N-SWNT-S is 
greatly enhanced and is much larger than the pumped 
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current / of the same sjistem when the superconduct- 
ing lead becomes normalEj. Roughly speaking, we have 
J^s ~ + Ea exp(i6') P where 6 is the phase difference 
between Ed and Ea- In the right inset of Fig. 2, we plot 
the ratio of jj^ as a function of pumping amplitude 
Vp at the resonant energy. At small pumping amplitude, 
the direct reflection and the multiple Andreev reflection 
are exactly in phase {9 = 0) and the complete,construc- 
tive interference gives the ratio around fouill3. As the 
pumping amplitude increases, the constructive interfer- 
ence effect in the N-SWNT-S system is suppressed be- 
cause 6* ^ 0. 

Now we consider a zigzag (10,0) SWNT with L = 56 
layers of carbon atoms (total 560 atoms). Two pump- 
ing driving forces are added on the tube layers from 
5th to 8th layer and from 49th to 52nd layer. By ad- 
justing Vq ~ 2.10065y, one resonant level is available 
at Ep = 0. In the calculation, we set (j> = 7r/2 and 
Tr^Tl ^ 0.0136eV". The results are plotted in Fig.3. 
We see that the pumped current gives similar behavior 
to that of armchair SWNT system. The pumped current 
is large near the resonant level and also has the double- 
peak structure around the resonant point. The ratio of 
jNS I jN jg ^Yso about four in the weak pumping regime 
(see the inset of Fig.3). Moreover there are several points 
worth mentioning. First, the pumped current can either 
be positive or negative even for the same SWNT but 
different energy levels or different phase differences (see 
Fig. 5 and Fig. 6). Comparing with armchair structure, 
the pumped current for zigzag structure reverses the di- 
rection. This is because the pumped current is not due 
to the external bias but cause by the pumping poten- 
tials. As a result, the pumped current is very sensitive 
to the system parameters. Second, the double-peak for 
the pumped current is asymmetric especially when the 
pumping amplitude is large (see the dashed hne of Fig.3). 
This is mainly due to the energy dependence of tbe. self- 
energy. Third, under the same system parameter£l, the 
pumped current of zigzag structure is much larger (at 
least ten times larger) than that of armchair structure. 
This means that in order to obtained large pumped cur- 
rent, one should use zigzag nanotube instead of armchair 
nanotube. This may be useful for experimental study of 
the carbon nanotube pump. 

Fig. 4 gives the pumped current as a function of pump- 
ing amplitude Vp at Ep = for armchair SWNT. Here 
the system parameters are the same as those in Fig. 2. We 
see that the pumped current increases quadratically at 
small pumping amplitude and then reaches linear regime 
for large pumping amplitude. In order to understand this 
figure, we also plot / = —1.06 x 10''!^^ in the same plot. 
We confirm that in the weak pumping regime the pumped 
current is proportional to the square of the pumping am- 
plitude, but in the strong pumping regime the pumped 
current is linearly proportional to Vp only. Similar con- 
clusion can be drawn from the zigzag nanotube (see the 
inset of Fig. 4). Fig. 5 and the inset of Fig. 5 present the 
pumped currents of armchair SWNT versus phase differ- 



ence at different pumping amplitudes in the weak pump- 
ing regime and the strong pumping regime, respectively. 
We see that the pumped current is anti-symmetric about 
(j) ~fl, just like the the result given by the experiments of 
Ref.cl although superconducting lead was not used there. 
In the weak pumping regime, the sinusoidal behavior is 
clearly seen. In the strong pumping regime, however, we 
see strong deviation from the sinusoidal behavior (see in- 
set of Fig. 5) and the maximum pumped current occurs 
near (f> = tt instead of 4>/2. 

Finally, we examine another N-SWNT-S quantum 
pump using a zigzag (10,0) tube with L = 92 layers (to- 
tal atoms 920). One gate is located from 10th to 28th 
layer and the other is located from 65th to 83rd layer. 
By adjusting Vq « 2.5936V, we obtain two double de- 
generated resonant levels at Ei = —3.5495 x lQ~^eV and 
E2 = 3.5495 X lO^^eT^. Hence, large Andreev reflections 
can occur near Ep — Ei and Ep = E2 with transmission 
coefficient equals to two (see the upper panel of Fig. 6). 
Here the Andreev reflection is due to different origin from 
that of Fig.l. In Fig.l, we have Eq — and Eq.([l^) gives 
us the resonance a.t E = with Ta ~ I- li Eq is nonzero 
the maximum Andreev reflection is less than one. In the 
strong tunneling regime which applies to our case, Fi 
and r2 are very small. To simplify the discussion, let us 
assume Fi = r2 << 1, then from Eq.(|l0|), we have 



^ [4(£;2-£;2)2 + rf + 4i;2r2] 

At resonance, Ta = r^/[r^ + AE^]. Hence when Eq is 
larger than Fi, the Andreev reflection quickly goes to 
zero. However, if the chemical potential of the supercon- 
ducting lead (^s=0 in our case) is right in the middle of 
two resonant levels {Ei and E2), i.e., fis = {Ei + £'2)/2, 
then electron coming from normal lead with incident en- 
ergy El tunnels into the structure through the resonant 
level El and Andreev reflected as a hole back to the quan- 
tum dot through the resonant level E2 with a Copper pair 
created in the superconducting lead, giving rise to the 
complete transmission. This is why in the upper panel 
of Fig. 6 we have two transmission peaks with Ta — 2 
due to the double degeneracy. The pumped current as 
a function of Fermi energy is plotted in Fig. 6. For com- 
parison, we also plot the transmission coefficient (long 
dashed line for normal structure and solid line for NS 
structure) in the upper panel of Fig. 6. The pumped cur- 
rent with = 1 X 10~^V^ (long dashed line in Fig. 6) 
for the N-SWNT-N system with the same system pa- 
rameters is also shown. Similar to Fig. 2 and Fig.3, the 
pumped current also shows strong resonant behavior. It 
has large value near energies where the Andreev reflection 
peaks occur, while it diminishes quickly away from the 
peaks. The amplitude of pumped current also increases 
as the pumping amplitude Vp increases. The striking fea- 
ture here is that the pumped current peaks have oppo- 
site sign for the two energy levels. That means the pump 
has property that the DC current can flow out of the 
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device from either electrodes by a slight change of elec- 
tron energy. From the lower inset of Fig. 6, we see that 
the pumped current clearly consists of two asymmetric 
peaks. Finally, we notice that large pumped current is 
generated for very small pumping amplitude (compare 
Fig. 3 with Fig. 6). This is because in Fig. 6 the thickness 
of potential barrier and hence the effective potential bar- 
rier height is much higher than that in Fig. 3. As shown in 
Ref. that the maximum pumped current can reach 

l/27r for extremely large barrier and in strong pumping 
regime. 

In summary, we have investigated the parametric 
pump of N-SWNT-S systems. By comparing with the 
parametric pump of N-SWNT-N system, we find that in 
the presence of superconducting lead, the pumped cur- 
rents is greatly enchanced due to the quantum interfer- 
ence of direct reflection and multiple Andreev reflection. 
In the weak pumping regime, the pumped current is pro- 
portional to the square of the pumping amplitudes but 
in the strong pumping regime, the dependence becomes 
linear. Hence large pumped current can be generated 
by increasing the pumping amplitude. When two level 
Andreev reflection occurs, the pumped current show re- 
markable parity effect so that the pumped current at one 
resonant level has opposite direction of that of the other 
resonant level. Our numerical results also show that the 
zigzag nanotube is a better candidate for the pumping 
device since it is more easily pumped than the armchair 
nanotube. In view of the S-SWNT-S structure studied in 
Ref. 
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it is conceivable that the N-SWNT-S quantum 
molecular pump can be realized experimentally and the 
results presented here be verified. 
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FIG. 1. Andreev reflection coefiicient Ta as a function of 
Fermi energy for N-SWNT-S system (solid hne) and transmis- 
sion coefficient T as a function of Fermi energy for N-SWNT-N 
device (long dashed line). The SWNT is au armchair (5,5) 
metallic tube. Inset: a schematic plot of the molecular device. 



FIG. 2. The pumped current / versus Fermi en- 
ergy for N-SWNT-S device at different pumping amplitudes: 
Vp = 0.0014V (sohd hno), Vp = 0.002V (dotted Unc) and 
Vp = 0.005F (dashed hne). Here long dashed line is the 
pumped current versus Fermi Energy for the correspond- 
ing N-SWNT-N device (when A = 0) at = 0.0014^. 
The left inset: j'^s ^^^.^^^ resonant point. The 

right inset: Andreev reflection coefficient Ta versus Ef with 
Vp = 0.005V at different pumping moments: t = 3n/4 (solid 
line), t = 0.95 x 37r/4 (dotted line) and t = 1.05 x 37r/4 (dashed 
line). 



FIG. 3. The pumped current /^"^ versus Fermi En- 
ergy for N-SWNT-S device at different pumping amplitudes: 
Vp = O.OOIV (solid hne), Vp = 0.002V (dotted line) and 
Vp = 0.005V (dashed line). The long dashed line is for the 
pumped current 7^ versus Fermi Energy of the correspond- 
ing N-SWNT-N device at Vp = O.OOIV. Here the SWNT is a 
zigzag (10,0) nanotube. The left inset: /'"^■^'//^'^ versus Vp at 
resonant point. 



FIG. 4. The pumped current /^"^ of the armchair SWNT 
as a function of pumping amplitude Vp at Ef = (solid 
line). Other system parameters are the same as those in Fig.2. 
Dotted line is plotted according to / = —1.06 x 10*1^^. Inset: 
The pumped current 7^'^ of the zigzag SWNT as a function 
of pumping amplitude Vp at Ef = (solid line). System 
parameters are the same as those in Fig. 3. Dotted line is 
plotted according to 7 = 8.3758 x 10^ Vp. 



FIG. 5. The pumped current 7'^'^ of the armchair SWNT 
as a function of phase difference (j> at Ef = 0. Main figure: 
the weak pumping regime with Vp = O.OOOIV (solid line) and 
Vp = 0. 0002V (dotted line). Inset: the strong pumping regime 
with Vp = 0.0014V (solid line), Vp = 0.002V (dotted line), 
and Vp = 0.005V (dashed line). Other system parameters are 
the same as those in Fig.2. 



FIG. 6. The pumped current 7 versus Fermi En- 
ergy for the N-SWNT-S device at different pumping ampli- 
tudes: Vp = 1 X 10~'^V (solid hne), Vp = 2 x lO^^V (dotted 
line). The long dashed line is the pumped current 7^ ver- 
sus Fermi Energy for the corresponding N-SWNT-N device 
at Vp = 1 X lO^'^V. The SWNT is a zigzag (10,0) nanotube 
with length L=92 layers. The upper inset: Andreev reflection 
coefficient Ta for N-SWNT-S device (solid line) and transmis- 
sion coefficient T for the corresponding N-SWNT-N device 
(long dashed line). The lower inset: the amplified figure of a 
pumped current peak. The other system parameters: (j> = 7v/2 
and Til = Pi, = 0.0136eV. 



6 




SNi luejjno pediund 




SNi luejjno pediund 




SNi luejjno pediund 




SNi luejjno pediund 




SNi luejjno pediund 



